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Abstract: We review the nuclear and condensed matter physics underlying the thermal 
and transport properties of the neutron star inner crust. These properties play a key role in 
interpreting transient phenomena such as thermal relaxation in accreting neutron stars, su- 
perbursts, and magnetar flares. We emphasize simplifications that occur at low temperature 
where the inner crust can be described in terms of electrons and collective excitations. The 
heat conductivity and heat capacity of the solid and superfluid phase of matter is discussed 
in detail and we emphasize its role in interpreting observations of neutron stars in soft X- 
ray transients. We highlight recent theoretical and observational results, and identify future 
work needed to better understand a host of transient phenomena in neutron stars. 
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1. Introduction 

Left undisturbed, neutron stars quietly fade away on a time-scale of about a million years 
as neutrinos from the core and, later, photons from the surface sap their thermal energy. 
The neutron star core, mysterious as it still remains, is expected to be a good conductor 
of heat, and with sufficient mass and heat capacity it acts as an excellent heat reservoir. 
In contrast, the crust is composed of matter that can be activated by external perturbations 
such as accreting matter, evolving magnetic fields and angular torques. Several observed 
transient phenomena such as X-ray bursts, superbursts, magnetar flares and glitches in the 
spin evolution are thought to originate in the crust. X-ray bursts occur at relatively shallow 
depth, at low density and high temperature, while superbursts, flares, glitches and heating 
during accretion are expected to occur deeper in the crust, in denser and more degenerate 
conditions. In these latter phenomena, thermal and transport properties differ from those 
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expected in ideal gases because of strong Coulomb and nuclear interactions, matter degen- 
eracy, and superfluidity in the inner crust. Interestingly, these properties play an important 
role in shaping temporal aspects of the observed X-ray light-curves. Recent observations 
of transient phenomena in accreting and magnetized neutron stars have renewed interest in 
developing realistic models of the neutron star crust. 

In this chapter we discuss the basic physics of the inner crust from the viewpoint of a 
low energy theory of electrons and phonons. For a more comprehensive survey we refer the 
reader to an excellent review article by Chamel and Haensel HJ. 

2. The composition of the inner crust 

The neutron star crust is composed of matter at sub-nuclear density (p < 10^^ g/cm^) and 
occupies the outermost 1-2 km region of the star. In the outer crust, matter is composed of 
fully pressure ionized nuclei immersed in a degenerate electron gas. With increasing den- 
sity, the electron Fermi momentum ^pe = (Stt^Z ^/)^/^, where Z and tii are the charge and 
number density of nuclei, increases rapidly and the reaction e~ + p ^ n + Ve drives nuclei 
to become increasingly neutron rich. Eventually, at p 4 — 7 x 10^^ g/cm^ the neutron 
fraction exceeds that which can be bound in the nucleus and some neutrons drip out and 
become delocalized from nuclei. The fraction of unbound neutrons increases rapidly with 
density and they become superfluid due to strong attractive interactions. The illustration in 
Fig. [T] shows typical variation in composition, density and electron chemical potential in the 
neutron star as a function of depth from the surface. Here it is assumed temperature is small 
compared to the energy differences between the ground state nucleus and other excited or 
metastable states. Under these conditions a unique ground state nucleus exists at each depth 
and is referred to as cold catalyzed matter. The region between neutron-drip, where the fa- 
vored ground state nucleus is ^^^Kr, and the core of the neutron star is called the inner curst 
and was first studied in pioneering work by Baym, Bethe and Pethick [2 J and by Negele 
& Vautherin [3|. In the shallower regions of the inner crust, spherical neutron-rich nu- 
clei are embedded in a neutron superfluid, but with increasing density non-spherical shapes 
(collectively referred to as "pasta" phases [4J) are likely to appear (see inset in Fig.[T]). 

The composition and structure of the inner crust have been extensively studied using 
mean field models based on the Skyrme interaction or its relativistic generalizations. These 
studies have elucidated how the structure and composition are determined by two key driv- 
ing forces: (i) the density dependence of the nuclear symmetry energy; and (ii) the gradient 
or surface energy in heterogenous asymmetric matter. Although the suite of models have 
explored a range of predictions for the symmetry and gradient energy, these ingredients 
remain poorly constrained by both theory and experiment. Hence the spatial extent and 
composition of the denser regions where pasta phases are expected, are active areas of re- 
search and is discussed in several contributions in this book. 

The composition in accreting neutron stars can be more complex since the crust here is 
formed almost entirely from accreted matter. Nuclear reactions that process this material 
have been studied in the literature and several relevant reaction pathways are now known 
[see Em], but the rates of some crucial pycno-nuclear fusion reactions at the temperatures 
attained in the crust remain very uncertain|7 |. Thus, several metastable nuclei can coexist 
at the same depth in the crust depending on its formation history. In the outer crust, the 
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Figure 1. Neutron star profile emphasizing the crust. The top x-axis shows the depth from 
the surface and left and right y-axes show the density and electron chemical potential, re- 
spectively. 

diversity in nuclear species is expected to be quite large in systems where nuclear burning 
and rapid-proton (rp) capture reactions during an X-ray bursts can produce nuclei with 
large Z and A . With increasing depth a sequence of electron capture, and neutron emission 
and absorption reactions tend to reduce both the diversity and the average Z and A of the 
nuclei encountered |[5l. These nuclei with lower Z and A are metastable because the fusion 
reactions needed to process them to the ground state are slow. 

The composition of matter is shown in the left panels of Fig. [2] where Z is the nuclear 
charge, A = A/^ + Z where N is the number of neutrons in the nucleus with single particle 
energy < 0, and Aceii = A/^+Z+ A^out is the total number of nucleons per unit cell (containing 
one nucleus) and A^out is the number of neutrons with single particle energy > in the cell. 
The panels show: (A) cold catalyzed matter where matter is in the ground state |3 1 ; (B) 
matter in accreting neutron stars where nuclear reactions including pycno-nuclear fusion 
process ^^Fe at p 10^ g cm~^ to higher density lO; and (C) matter in accreting neutron 
stars where pycno-nuclear fusion is assumed to not occur for Z > 4; are shown. In models 
(B) and (C) we use the composition of |3| for catalyzed matter for p > 2 x 10^^ g/cm^. In 
these plots we assume Z changes abruptly while A and Aceii vary smoothly, and the vertical 
dashed lines show the point at which neutron drip occurs. 

Unless stated otherwise, in this article we use natural units where the speed of light 
c = 1, Planck's constant /i = 1 and the Boltzmann constant = 1. Typically, we will 
quote energy in MeV and length in fermi and velocities in units of c. To convert to cgs 
units, quantities of interests such as the thermal conductivity and the specific heat per unit 
volume, we note that MeV^ = 1.07 x lO^^f-^l and MeV^ = 1.81 x 10^^[^^]. 

' Lcm s KJ Lcm-^ K-i 
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Figure 2. Left Panels: Crust chemistry. The variation of Z, A and Aceii with density. The 
vertical dotted line show the neutron drip point Pdrip. Right panels show characteristic 
temperature scales discussed in ^3]J Tm is the classical ion melting temperature, for a critical 
r 2^ 200, Tp is the ion plasma temperature, Tum is the crystal umklapp temperature, and 
is the ion de Broglie temperature. 



3. Low energy excitations and thermal properties 

To describe the thermal properties we will dwell on a key simplification of dynamics at low 
temperature where only electrons and collective excitations or phonons are relevant. As we 
shall see below the details of the complex nuclear ground state and the lattice structure are 
encoded in a few low energy constants which determine the propagation and interaction of 
phonons. Three of these modes are phonons of the lattice, and can be thought of a coherent 
motion of the clusters of protons that are located at the center of the neutron-rich nuclei. 
The forth mode is associated with neutron density oscillations and is called the superfluid 
phonon mode when neutrons are superfluid. The motion of proton clusters and neutrons are 
not independent, and a fraction of the neutron fluid is entrained in the motion of the proton 
clusters. The number of neutrons entrained by each proton cluster is in general different 
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from N and as we shall discuss below it is a key microscopic input for the low energy 
theory of phonons. 

To begin we shall briefly discuss the separation between the low and high energy scales 
for each of the three components in the inner crust shown in Fig. [3] 



3.1. Ions 

Since the electron screening momentum ^tfc = \/4e^/7Z ^Fe ^ {2a)~^ where 2a is inter-ion 
distance, ions crystallize at T ^ — Z^e^/{rM ci), where Tm is the melting temperature 
of the Coulomb solid. For a one component Coulomb system Fm = 173 but can be as large 
as 230 for a mixture a elements. Written in terms of the chemical potential /i^ = ^pe for 
relativistic electrons, 



which implies that the inner crust where /i^ ^ 30 — 100 MeV is a solid for T < 10^ K. For 
ions, the intrinsic energy scale is the plasma frequency 00^ = {A%e^ni 7? jMf)^l^ where 
/!/ = 3/471 is the number density of ions, and Z and M\ are the charge and effective mass 
of the ions, respectively. In terms of /i^. 



CO, = = ^ 0.3 f ^^^)'^VeV, (2) 

^ ^ y 371 M; Va* V30MeV/ ' ^ ^ 

where A* jrun is the effective number of nucleons that move (entrained) with each 
nucleus. 

Below the Debye temperature Td ^ 0.457, one longitudinal and two transverse phonons 
of the lattice capture all of the relevant ion dynamics. Uncorrelated single particle motion 
for ions does not emerge until T ^ Tp. The variations of Tm and Tp with density and 
composition are shown in Fig. [2] Generally ,we can expect that nuclear excitations are 
suppressed by shell gaps and pairing correlations. The evolution of these effects in the 
inner crust is yet to be understood. We denote this as 7a ^ 1 MeV - a typical nuclear 
physics scale. It is possible that surface modes and low-lying collective excitations in the 
large nuclear structures encountered in the pasta-phases could have smaller energy. 



3.2. Electrons 

In terrestrial solids, the electron excitation spectrum is complex as band structure and 
Cooper pairing play a role at low temperature manifesting in insulating, metallic and super- 
conducting properties. In contrast, electronic structure in the crust is simpler because the 
electron Fermi momentum ^pe is larger than all other relevant scales depicted in Fig. [3] 

Electrons are relativistic with v/7 2^ 1 and the ratio Ve-i/A'g — ^^^^ <^ I where 
Ve-i = Ze^/a is the interaction energy and He = (971 Z/4)i/V a is the electron chemical 
potential. Consequently the Fermi surface is nearly spherical with minor distortions at in- 
tersections with Brillouin zone boundaries where coherent Bragg scattering leads to a band 
gap[8|. These small patches on the Fermi surface become important only at low temper- 
ature when considering phonon mediated transitions between electrons. As we discuss in 
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Figure 3. The regimes where collective excitations dominate over single particle excitation 
for ions and neutrons. Scales of relevance to electron dynamics are also shown. 

^4} the gap affects the Umklapp process for T £ Tum = Vt hu where hu is the band gap 
and Vt is the velocity of transverse phonons |9|. For nearly free electrons hu = Vkpe^ where 
Vkp^ ^ AnZe^ nj/k^^^ (4e^/3K) kpe is the Fourier component of the lattice potential at scale 

Although electrons in the inner crust are as degenerate as terrestrial superconductors 
with T/Tp ^ 10~^ — 10""^, where Tp = jUe, the critical temperature is negligibly small be- 
cause here electrons are relativistic. They move too quickly to adequately experience the 
attraction due to retardation effects in the electron-phonon potential, and consequently the 
critical temperature T/ 2^ co^ exp {—vp/e^) <C COp is negligibly small |[TOll . Thus, the degen- 
erate Fermi gas model provides an excellent description of electronic properties for T <Tf. 
In this regime, the density of states A^g(O) = /i^/ji^ is large and this greatly enhances their 
contribution to thermal and transport properties at low temperature. 

3.3. Neutrons 

Due to strong attractive interactions, neutrons in the inner crust form Cooper pairs and be- 
come superfluid. The gap in the single particle spectrum is denoted by increases from 
zero at neutron drip to a maximum value 1 MeV at a density p 2^ 10^^ g/cm^ and de- 
creases therafter. The number of thermally excited neutron quasi-particles is exponentially 
suppressed when T <C A^ and their contribution to thermal and transport properties is typ- 
ically negligible. However, depending on the variation of the gap with density, a sizable 
fraction of the inner crust close to neutron drip and the vicinity of the crust-core interface 
can be normal in accreting neutron stars. 
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Figure 4. A sample of theoretical prediction for the neutron ^Sq superfluidity critical tem- 
perature Tc. 



In Fig.|4]model predictions for the critical temperature Tc = A^/1.76 are shown where 
curves labelled "BCS" and "GMB" show the analytical results in the weak coupling valid 
in the limit lakp] ^ 1. In the Bardeen Cooper and Schrieffer (BCS) approximation Abcs = 
(8/^^) exp(7i/2a kf)EF, with a scattering length a = —18.5, fm. Corrections due to medium 
polarization which appear at the same order reduce the gap to Agmb = 1/ (4^) Abcs from 
flu . Curves labelled "Al" and "A2" are examples of slowly growing Tc at low kp, from fl^ 
and 1 13], respectively. Curves "Bl" and "B2" mimic behavior predicted by strong coupling 
QMC calculations from [14] and [TSl where the gap increases rapidly with density. In 
models labelled "Al" and "Bl" where gaps vanish at p 2^ 10^"^ g/cm^. For more details on 
the density and model dependence of the gap we refer the reader to the chapter by Gezerlis 
and CarlsonCH in this book. 

In the region where T <Tc collective excitations of the neutron fluid called superfluid 
phonons, with a dispersion relation co = V(|) ^, are the relevant low energy degrees of free- 
dom. This mode corresponds to fluctuations of the phase of the superfluid condensate (and 
can be related to density fluctuations) and is the Goldstone mode associated with the spon- 
taneous breaking of the global U{\) symmetry in superfluid ground state (the Hamiltonian 
is invariant under arbitrary phase rotations of the fermion fields, but in the superfluid ground 
state is preserved only by discrete rotations of 7i/2). 



3.4. Specific heat 

The electron contribution the specific heat (hereafter Cv will represent the specific heat per 
unit volume) is given by 

Cl = \nlT, (3) 

at low temperature. Band structure affects only negligible as only small regions of the 
Fermi surface are affected. At low-temperature when T <^Tp electrons dominate, but as 
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we discuss below the phonon contribution can become important in accreting neutron stars 
where T 10^ — 10^ K. For T the contribution from lattice phonons (Iph) is given by 

cr=t(K^'^). (4) 



15 Vv^ / 

and V/ and Vt are velocities of the longitudinal, and transverse lattice phonons, respectively. 
In a model where the strong interaction between the neutron superfluid and the ion lattice is 
ignored it is simple to calculate these velocities. The speed of longitudinal lattice vibrations 
is approximated as v/ = y^Kio^-^/p where iiTion-e = p(9(/^ion + /9p) is the bulk-modulus 
of the electron-ion system and the ion mass density p = Arun nj where A is the number of 
bound nucleons in the ion. Since ^ ^ion? we can write 

which is usually referred to as the Bohm-Staver sound speed. The velocity of the transverse 
lattice mode is related to ju, the shear modulus of the lattice, and is given by 

V, = J^ = a^, (6) 
VP 

where = (671^^/)^/^ is the ion Debye momentum, and the constant a c:^ 0.4 is obtained 
by numerical calculations of Coulomb crystals fTfi . Further, since 




^TFe 




» 1 (7) 



we have v/ » Vt and the contribution from longitudinal modes to Cv in Eq.|4]in negligible. 
Thus the lattice contribution can be written in the familiar form 

where To = (3/2) ^/^v^ — 0.45 is the Debye temperature of the ion lattice. This low 
temperature form of the specific heat provides an excellent approximation in Coulomb 
solids up to r < Tp/50 but fails when T > T^/IO 1 18 1. 

To calculate the neutron contribution to Cv we first note that there are two distinct 
regimes. In the normal phase when T >Tc the neutron contribution is large and is given by 

C~ = ^m„^F„r iT>T,) (9) 

This normal contribution can become important in the vicinity of neutron drip where T >Tc, 
and at the crust-core boundary. In the superfluid phase when T <^Tc the neutron single 
particle excitations are strongly suppressed and 

-^neutron _ ^ t rr. / 



c~"~3"tn^F„rexp(^y^j (r«r,). (lO) 

which is usually negligible. The four models for the gap in Fig. |4] allow us to explore the 
effect of pairing on the neutron specific heat. In models Al and A2 we have a thick shell 
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of normal neutrons above the drip point, while models Al and Bl predict a thick layer of 
normal neutrons at the highest densities. Modifications to this simple picture of pairing in 
uniform neutron matter due to the presence of the nuclei are discussed in this book in the 
chapter by N. Sandulescu & J. Margueron |[T9ll . Further, we briefly note that like in the 
case of electrons, coherent Bragg scattering of neutrons by the lattice lead to band structure 
effects that modify the shape of the Fermi surface, still Eq.|9]is an excellent approximation 
to Cv in normal phase for reasons described in [20]. 

Elsewhere in the crust where T <^Tc the relevant neutron contribution is from superfluid 
phonons, i.e., collective instead of single particle excitations, and is given by 



C; 



(Sph 



15^ 



(11) 



where 



with/^2 



see ^331 



(12) 



is the superfluid phonon velocity, rif, jUn and rrin are the number density, chemical potential 
and mass of the free neutrons, respectively. For weakly coupled systems = vf/\/3 where 
vp is the Fermi velocity. In most of the inner crust » Vt (see Fig. [6]) and hence their 
contribution to the heat capacity is negligible except perhaps in a sliver where ^ and 
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Figure 5. Specific heat of ions, electrons, and for neutrons with (labelled ^sf) and without 
the effects of the superfluid gap (labelled n^) are shown for four representative tempera- 
tures. 

The specific heat due to these components is shown in Fig. [5] The ion contribution for 
T ^O.lTp varies as and is to very good approximation given by Cy^, while electron 
contribution is linear in T and dominates at low temperature. As mentioned earlier, the 
neutron contribution is sensitive to the variation of the ^Sq gap. To illustrate this we show 
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predictions representative of two distinct possibilities corresponding to the behavior akin to 
models A's and B's of Fig. |4j In models labelled Al and A2 the neutron contribution is rele- 
vant over an extended region near neutron drip, while in models Al and B 1 the contribution 
is relevant in the region close to the crust-core boundary. The contribution of the superfluid 
mode is negligible and is not shown. The ion contribution in the phonon regime T < To 
changes quite significantly with the effective mass of the nucleus since Cy^ oc Mj^^'^: its 
increase from a = to 30% and 60%, where a is the fraction of unbound neutrons entrained 
by each nucleus, is illustrated in Fig. |5j We now turn to briefly describe the low energy 
theory where effects due to entrainment are taken into account consistently. 



3.5. Entrainment and mixing 

In the preceding discussion the interaction between the neutron superfluid and the lattice 
was ignored. Because neutrons and protons interact strongly this is difficult to justify and 
recent work in ll2Tll describes a first attempt at including these interactions. Using the tech- 
niques of low energy effective field theory, a theory of lattice and superfluid phonons was 
formulated where the relevant fields (]) and correspond to the phase of the neutron field 
and the displacement vector of the proton clusters from their equilibrium positions, respec- 
tively. The symmetries associated with translation and number conservation require that the 
low energy theory be invariant under the transformation (r, t) (j^ _^ a^^^ -^ 

and (l)(r,/) (l)(r,0 + 9 where a^^^-^ and 9 are constant shifts. Since the theory must re- 
spect these symmetries, it implies that the Lagrangian can contain only spatial and temporal 
gradients of these fields. This gradient expansion enables us to organize the calculation of 
any low energy process in a well controlled expansion in powers p/A where p is the mo- 
mentum (involved in the process) and A is an intrinsic high-momentum scale associated 
with the ground state. Assuming that the ground state must posses a cubic symmetry, the 
quadratic part of the effective Lagrangian density 
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ab 



where higher order terms involve higher powers of the gradients of these fields, and ^' 

If the superfluid and the solid are uncoupled, the low energy coefficients (LECs) appear- 
ing above, such as p/,/i/, ilT/ are the mass density, the shear modulus, and the compressibility 



of the solid, respectively, and and of the superfluid are as defined earlier in Eq. 12 
They determine the velocities 



v/ = W , V, = , and Z^-^, (14) 

VP/ VP/ Y fi 

respectively, where Uf now corresponds to the number density of free neutrons and p/ = 
Anin nj the mass density of the nuclei. In the presence of strong coupling between the solid 
and superfluid these coefficients are modified. As the proton clusters move, they drag along 
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neutrons from the fluid and this will modify the mass density involved in lattice fluctuations. 
Correspondingly, the coefficient p/ appearing in Eq.[T3]and the lattice velocities will differs 
from the usual mass density of the pure lattice component due interactions that entrain 
the superfluid ||22l[23l. Naively, one may associate the number of neutrons entrained by 
each proton cluster to correspond to the number of bound neutrons with single-particle 
energy less than zero. However, this typically underestimates the mass-density associated 
with lattice motion. Calculations indicate that a large fraction of neutrons with positive 
single particle energies are entrained due to coherent Bragg scattering of neutrons from 
the lattice [24 J . Denoting this number density of entrained neutrons as nt we can write 
Pi = {nb + np)mn where rip is the average proton number density in the unit cell. 

It is convenient to express njy in terms of the neutron effective mass nf calculated from 
band structure studies outlined in 
j% 



Ace,i(l-^) +a'^-Z={A*-Z), (15) 



where Aceii — p/mn is the total number of nucleons in the unit cell, A is the number of 
nucleons bound in the nucleus and A* is the effective number of nucleons that move with 
the nucleus: 

TTl 

A*=A + a(Aceii-A) with a=l- — . (16) 

K 

Representative values of Aceii, A, Z in the inner crust as function of total mass density p 
are shown Fig. [2]. Calculations reported in [241 indicate that m* is typically in the range of 
(3 — 6) Mn but can be as large as 15 is some regions where coherent Bragg scattering 



is most efficient. From Eq. [T5] we see that when m* » a large fraction of neutrons are 
entrained by nuclei. 

Current conservation implies that the number of neutrons that move freely as part of the 
superfluid is correspondingly reduced and this is denoted by ^/ = — ^z?) where rin is the 
total neutron density [25 1. This repartition of neutrons between the lattice and the superfluid 



will modify the LECs appearing in Eq. 13 Now p/ = (nt + np)mn and Uf = (un — nt) 



in Eq. [T4j The thermodynamic derivates needed to define the lattice compressibility K 
the shear modulus ju also contain the effects of the neutron-proton interactions as they are 
accounted for in the calculation of the ground state properties. The specific thermodynamic 
derivatives that define the LECs of the mixed system are explicitly given in f2r|. 

Interactions between neutrons and protons leads to mixing between the modes of the 
lattice and the superfluid. The extent of this mixing is controlled by the dimensionless 
coefficient 



1 / . 

fib-rip^ , (17) 



^mix = 7 — I rib - np^"" 



where the first term in the parenthesis arises due to the interaction between neutron and pro- 
ton velocities (entrainment) and second term is due to the interaction between neutron and 
proton densities [21]. Estimates indicate dfin/drip Krib/rip and mixing is mostly determined 
by the term proportional to ni, ||25l[2Tl and approximately we can write 
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Figure 6. Velocities of phonons in the inner crust for two chemical compositions from Fig 
[2] and three values of the nuclei effective mass A*. 

Mixing implies that the longitudinal eigenmodes are superpositions of the longitudinal lat- 
tice and superfluid phonons. The velocity of these eigenmodes is given by 



VI ,2 



l±l/l-^| (19) 



where X = g^-^ + v^ + v^ ^^d v/ and are defined in Eq. 14 The velocity of the eigenmodes 
for the crustal compositions of catalyzed and accreted matter shown in panels (A) and (B) of 
Fig.[2]are plotted in Fig.[6j The dashed curves show results for v/ and without mixing and 
they cross at p ^ 10^^ g/cm^. In this resonance region mixing is large and level repulsion 
can be significant. Away from resonance, the eigenmodes contain only small admixtures: 
below p 2^ 10^^ g/cm^ the mode labelled V2 is predominantly the superfluid mode and above 
it is predominantly the lattice mode. In these calculations we have neglected the second 
contribution in Eq.[T7]to g^ix and the value of nt was chosen somewhat arbitrarily to reflect 
the range of m* predicted in [24]. The panels show results for three values of g^ix chosen to 
reflect different fractions a = 0, 30% and 60% of unbound neutrons in the cell entrained by 
each nucleus. Transverse modes are unaffected by mixing at leading order but are affected 
by entrainment. Its variation in the crust for different values of a is also shown in Fig. |6] 
Despite strong mixing Vt <C vi or V2, and transverse modes will continue to be dominate the 
specific heat. 
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Figure 7. Feynman diagrams indicating the various scattering and dissipative processes 
involving electrons, lattice phonons and superfluid phonons. 



4. Transport Properties 

The electron and phonon thermal conductivity can be written as K = Cv v ?i/3 where Cv is 
their specific heat, v is their velocity, and X is the transport mean free path. Using Eqs.|3]& 
|4]the electron and phonon conductivities are 

1 o 271^ Q 

Ke = -Jill T Xe , Kph, = ^ (20) 

where electrons are relativistic (v = 1) with mean-free path Xe, and the phonon contribution 
is for each phonon type with velocity v/ and mean free path ?iph.. Since jLie ^ T, electrons 
dominate at low temperature but phonon contributions can become relevant at high tem- 
perature when ^ {/Ue/T)^ v? Xe or when the magnetic field is large enough to restrict 
electron motion I26,i27l. Phonon velocity was discussed in ^3]J we now turn to discuss 
scattering and absorption processes that determine their mean free path. Feynman diagrams 
for relevant interactions are illustrated in Fig. [7] and in the following we briefly discuss the 
most important of these processes in the inner crust. 



4.1. Electron-phonon processes 

In its general form, the electron mean free path relevant for the thermal conductivity due to 
electron-ion scattering is given by 

22^4 ^2kfe 

where 



K^ = -r^ dkk^\V{k)\^ JcoJ(p(o)S(co,A:)gK(A:,|3co) (21) 
4ji//| Jo J-oo 



and the dynamical structure factor 5(co,^) embodies all relevant dynamics of the strongly 
coupled system of ions ll28]| . Here, CO, k are the energy and momentum transfer. The function 
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Figure 8. Feynman diagram for the Umklapp process. The double dashed line represents 
recoil-free momentum transfer Q^k — qioihc lattice, and \q\ Kqo lies in the first Brillouin 
zone. 

V{k) = Fz{k) I (k^ + ^TFe) characterizes the screened electron-ion interaction in momentum 
space where ^^p^ = Ae^]^p^l% and Fz(^) is the charge form factor of the nucleus. 

Pauli blocking restricts CO ^ T <C /i^, and when 5(co,^) contains most of its strength in 
the region CO ^ 37 the conductivity can be expressed in terms of the static structure function 
S{]i) = / Jco S(k^ co). However, 5(co, ^) has strength at co ^ co^ and \e cannot be calculated in 
terms of when T <Tp. Here, the frequency dependence of the dynamic structure factor 
is needed but this is generally difficult to calculate in strongly coupled quantum systems. 
Fortunately, when T < To phonons are the only relevant degrees of freedom and electron 
scattering is dominated by the emission or absorption of phonons |29|. In this case, 5(co,^) 
is simpler and is characterized by discrete peaks at CO = associated with the excitation of 
phonons with velocity v. 

In the low-energy theory, the interaction between electron and phonons is described by 
the Lagrangian density 

"^-^^-i^jl"^^^'^' -here /eph = (23) 

is related to electron-phonon coupling constant [305, is the electron field and is the 
ion displacement (phonon) field discussed in ^3]] This form of the interaction applies to 
normal processes, where the momentum transfer k < qo and displacements correspond to 
excitation of longitudinal phonons. However, since kfe/qo = > 1 large angle 

electron scattering with k > qo possible. This Umklapp process is depicted in Fig. [8] 
where the electron simultaneously Bragg scatter off the lattice and excite a phonon. Elastic 
Bragg scattering (without phonon emission) however does not contrbute because electrons 
are eigenstates of the lattice potential. Further, unlike normal processes where only longi- 
tudinal modes are involved, Umklapp scattering is dominated by the emission or absorption 
of transverse phonons |[28ll3T1l . 

The dynamic structure factor for single-phonon emission and absorption including 
Umklapp shown in Fig. [8] is given by 

^^^''^i-M^LL 2(0 Vl-exp(-P(o)+exp(-|3co)-ly'^ ^ ^' ' ^^^^ 

where the first and second terms in parenthesis represent phonon emission and absorption, 
respectively |28 1 . The phonon momentum is restricted to the first Brillouin zone q < qo. 
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and sums are over all reciprocal lattice vectors or lattice momenta Q, and the longitudinal 



and transverse phonon states with polarization vector 8^ and velocity v/. Using Eq. 24 and 



the delta functions to perform the integration over k and CO, the electron mean free paths in 



Eq. 21 can be written as 



1 271^6^ CO^ 



^!^«(r,v,-), where, (25) 

(26) 



Q Jo (27i)3 



(exp(Pv/^) - 1)(1 -exp-(Pv,-^)) 



and ^ = q + Q. To unravel the dependence on the temperature and the phonon velocity we 
examine two limitings forms of the function ^^^\T,Vi). First, when 2kfe ^ ^d, the domi- 
nant contribution comes from the Umklapp and we can set 5^ = 2 in evaluating (r, v/). 
In this case, from the RHS of Eq. [26]it is easy to deduce that 

lim!^W(r,v,-)oc— . (27) 
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to find that 



In the opposite limit, when only the normal process involving longitudinal lattice modes 
contribute we can set 2 = in the RHS of Eq. 

Q=0 ^ ' ' vj 



lim!^»(r,v/)oc^. (28) 



At very low temperature, the band gap in the electron spectrum suppress Umklapp 
processes. As mentioned in ^3] coherent Bragg scattering by the lattice will distort the 
electron Fermi surface for momenta that can coincide with the reciprocal lattice vectors 
Q. Here, the spectrum will differ due to a band gap 5u — {4e^ /37z) kj^Q. Although distorted 
patches on the Fermi surface occupy only a small fraction of the total area, these regions are 
important for Umklapp transitions. To understand this suppression consider the case when 
the phonon momentum ^ ^ 0. In this limit, large angle electron Umklapp scattering with 
^ 2^ 2 can only involve electrons on these patches. However, at low temperature the gap 
will suppress such transitions unless the phonon momentum q>5k where 8^ ^ 8f//vFe can 
"steer" electrons away from these patches. For transverse thermal phonons qc:±3T /vt and 
the condition on the phonon momentum implies that Umklapp occurs for T > Tum where 

From the preceding discussions we can conclude that for T > Tum the mean free path 
ocv^/T^ since <C v/. For T <C T^m where only normal processes involving longitudinal 
phonons are allowed we expect Xf^ ocyf /T^. However, the normal electron-phonon process 
is too weak to compete with two other sources of electron scattering that we now discuss. 



4.2. Electron-impurity scattering 

As we noted in ^2] in accreting neutron stars nuclear reactions that process accreted mate- 
rial can produce a mix of metastable nuclei. The evolution of nuclei in the outer crust has 
been studied in ||5l where it was found that electron capture induced neutron emission reac- 
tions populate a very diverse mix of nuclei with a large dispersion in Z and A. Although it 
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is reasonable to expect that this dispersion will significantly decrease in the inner crust due 
to pycno-nuclear reactions and the abundant supply of neutrons, reaction pathways in the 
inner crust remain poorly understood. It is generally assumed that at each depth a specific 
nucleus with large Z and the highest abundance will crystallize and the remaining mix of 
nuclei can be treated as impurities in the solid. The impurity parameter 

amp = ^ I^^/(Z,■-(Z))^ (29) 

'^ion / 

is a good measure of the dispersion in the nuclear charge. For moderate Qimp ^ 1 an ordered 
lattice is likely with scattered impurities. If the impurities cannot diffuse easily their spatial 
distribution will be uncorrelated, and electron scattering off them can become significant. 
The scattering mean free path in this case is given by 

5, imp ^Fe a -1 ^TljZ) i 

471.4 j:.^.(Z,._(Z))2^ 4e^Q,mpk,, ' ^ ^ 

where A 2^ 1/2 (In (tc/.^) — 2) is the Coulomb logarithm, and we have used charge neutral- 
ity which requires (Z) ^lon = = kl^/37Z^ in arriving at the second equality. 

4.3. Electron-electron scattering 

Typically electron-electron scattering is weak but it can become important when electron- 
ion scattering is suppressed at T < Tum. Scattering between relativistic electrons is dom- 
inated by the current-current interaction which unlike the Coulomb interaction between 
charges, this interaction is unscreened in the static limit. The corresponding mean free path 
was calculated including the effects of dynamical screening (or Landau damping) in |[32ll . 
For the case of relativistic and degenerate electrons 

^^-^"6^[3] ^^^^ 

and it is remarkable that it is independent of density. The corresponding conductivity Ke_e — 
21 /i^ is also interesting as it is independent of temperature. Consequently, electron-electron 
process can become important at T < Tum when electron-phonon Umklapp scattering is 
suppressed. However, in practice for T > 10^ K they are only relevant in a small region 
close to the crust-core boundary if gimp <C 1. 



4.4. Electron conduction 

Numerical calculations of the electron conductivity with several refinements that include 
the role of multi-phonon excitations, Debye-Waller corrections and the nuclear form fac- 
tors have been calculated and tabulated by the neutron star research group at the loffe in- 
stitute in St. Petersburg (http://www.ioffe.rssi.ru/astro/conduct/I). Since our focus here is 
to emphasize the qualitative aspects at low temperature we do not review these important 
refinements. The results obtained (using the fits to the tabulated results) are shown in Fig. [9] 
and qualitative features can be generally understood in terms of our preceding discussion. 
Four panels with increasing T in Figj9] clearly demonstrates: (i) the rapid decrease in ther- 
mal conductivity for the case gimp = as T becomes larger than Tum and (ii) the importance 
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Figure 9. Electron thermal conductivity vs density at four different temperatures. Scat- 
tering processes e-ion, e-e, and e-imp with 6 values of Qimp = 0, 1, 2, 3, 4, and 5 (as 
indicated in the left panel) are included. 



of impurity scattering in the inner crust for T < 10^ K and for gimp ^ 1. Both of these trends 
are easily understood in terms of the preceding discussion of various scattering mechanisms 
and their temperature dependencies. As we discuss in ^5l| 2imp will play an important role 
in interpreting observations in accreting neutron stars when the inner crust is cold with 
T < T^m- 



4.5. Phonon conduction 



Phonon heat conduction can become relevant when T ^ 10 K when the phonon heat capac- 
ity becomes comparable to that of electrons, or when the electron contribution is suppressed 
either due to large Qimp or magnetic fields. Its importnace depends on the phonon mean free 
path being large enough to compensate for their smaller velocity. Phonon scattering pro- 
cesses have been discussed in Refs. |[26l [27l and we will briefly review them here. As in 
terrestrial metals |29|, electrons in the inner crust are efficient at damping lattice phonons. 
The phonon-electron process is shown in Fig [7] (2a) which depicts a phonon decay produc- 
ing an electron-hole excitation. This, Landau damping, dominates over phonon-impurity 
and phonon-phonon processes for the temperature realized in the crust l26|. 

The electron-phonon process discussed in §4.1. and the phonon-electron process we 
discuss here are essentially similar. Only here it acts to bring into equilibrium the phonon 
distribution function that carries the net thermal current relative to the electron gas. Since 
transverse modes dominate the heat capacity their contribution to thermal conduction is rel- 
evant and longitudinal modes can be neglected. For T > Tum, Umklapp processes dominate 
and transverse phonons are absorbed and emitted by large angle electron scattering on the 
Fermi surface. The mean free path for these processes was estimated by Chugunov and 
Haensel in ll26l . For simplicity, neglecting corrections due to the Debye- Waller factor, we 
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Figure 10. Electron and phonon thermal conductivity in the presence of large magnetic 
fields from |[27ll . 



can rewrite their estimate as 

Try 1 F{Tp/T) _ 283 a F{Tp/T) 



^Iph-e 

where 
F{Tp/T) 



V^Fe Aph-e vZ2/3 Aph-e 



= 0.014 + 



0.03 



exp(V5r) + l 



A 



ph— e 



2 V 4 



(32) 



(33) 



where v = (Op/3qj) Vt is average phonon velocity and Y= ^o/fe = (2/Z)^/^. For typical 
values of Z 2^ 20 — 40 and v ^ 10~^ the mean free path ^ 100 a and the corresponding 
thermal conductivity is negligible |[26l. Below Tum electron Umklapp is suppressed and 
phonon-electron process cease to operate for transverse phonons. Although their mean free 
path can be large, numerical calculations show electrons continue to dominate. Normal 
phonon-electron process will continue to operate for longitudinal phonons at T < Tum and 
their mean free path 



^Iph- 



71 



6C[3] 



1 

f 



(34) 



was estimated in ll27l . 

The superfluid can also transport heat. In terrestrial superfluids such a liquid helium 
extremely efficient heat flow occurs through a non-diffusive process called internal con- 
vection. This process which is also responsible for the propagation of second sound in 
superfluid helium is described as the counter flow of normal and superfluid components 
in the Landau two-fluid model |[33l . In neutron stars the situation was found to be quite 
different due to strong dissipation of the normal phonon component [27 J . Here, heat flow is 
diffusive as superfluid phonons are effectively scattered by their indirect coupling to elec- 
trons. This coupling is induced by the interaction between superfluid and lattice phonons 
described in |3.5.| and is shown in diagram (3a) of Fig. [7] The superfluid phonon mean free 
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path 



_ J_ l + (l-a^)^ (aniph)^ 



=»mix 



where a^vi/v^, Xiph = ?^iph-e/v/ and ?iiph-e is the mean free path of the longitudinal lattice 
phonon. Eq. [35] is valid in regions away from the resonance region where a ^ 1. Near 
neutron drip a » 1 and Eq.[T7]can be approximated as 

^ v] nf{nt^np) _v] (Aeeii-A*)A* 

Asph-e - 3 2 ^Iph-e " 3 - 7^2 ^IP^-e ' V-^O; 

to deduce that ?isph-e ^ ^iph-e- Superfluid heat conduction can become relevant here, espe- 
cially in the presence of large magnetic fields that suppress electron conduction transverse 
to the magnetic field. This is illustrated in Fig. [TOj With increasing density superfluid heat 
conduction becomes negligible because ?isph-e ^ ^iph-e in regions where v/ ^ and when 
V(|) ^ V/ the contribution remain small because Cv^^^^^ ©c 1/v? is negligible. 



5. Observable manifestations 

Possible observational probes of the thermal and transport properties of the inner crust re- 
quire astrophysical settings in which the crust evolves rapidly on observationally accessible 
timescales. A first case is simply the early cooling of a new-born neutron star, which relaxes 
from the initial hot proto-neutron star stage. The second family of scenarios involves some 
process which deposit a large amount of heat in the crust and observation of the relaxation 
of the crust can provide invaluable information. Such heating is expected to occur in a neu- 
tron star undergoing accretion from its companion in a binary system and, in a completely 
different setting, in magnetar giant flares where catastrophic magnetic realignment deposits 
a large amount of energy. In both cases thermal relaxation of a heated crust have been 
observed. There is also an intriguing possibility that the early evolution of a "neo-neutron 
star" born in the aftermath of a supernovae associated with GRBs has been observed in the 
late X-ray emission [i34i . 



5.1. Thermal evolution equations and crust relaxation 

The equations controlling the time evolution of the neutron star temperature T are 

dT 1 3L dT . 

Cv^ = Qh-Qv ^ and F = -k^ with L = 47irF (37) 

ot 47rr^ or or 

where the first equation simply express energy conservation and the second describes heat 

transport!^ F and L are the diffusive flux and luminosity, resp., within the star, Cy and K 

the specific heat and thermal conductivity discussed above, 2v is the neutrino emissivity 

and Qh the heating rate (both typically expressed in units of erg cm~^ s~^). Together these 

two equations give us the heat equation 

Cv^ = K-y + ^^-^— +e/,-ev (38) 

dt or^ r^ or or 



^ These are the Newtonian version of the full GR equations, see, e.g., 1351 . 
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from which a thermal time-scale can be estimated as 

Tth - ^ (39) 

where / is a typical length-scale for the T variation. So, the observation Xth would constrain 
Cy and K but only through the combination in Eq. [39} and modulo an assumed value of /. 
More information, however, can be obtained in special circumstances which we describe 
below. 

5.2. Crust Relaxation in Quasi-Persistent Accreting Neutron Stars 

Among neutron stars in accreting binary systems, a very special class called Soft X-Ray 
Transients (SXRTs) are formed by systems where accretion is transient. Among the 3 dozen 
known SXRTs, four systems have the peculiarity that the observed accretion outburst(s) 
lasted several years (while most SXRTs show accretion outbursts lasting a few weeks, sep- 
arated by months/years/decades of quiescence). These four systems have been dubbed as 
"Quasi-Persistent Sources" and they are excellent systems to study thermal relaxation. 

In these sources nuclear reactions in the crust deposit heat deep in the crust. The ac- 
creted H/He can be processed to large A 2^ 100 by the rapid proton (rp) capture process f36l 
and in the event of a superburst this mix is reprocessed iron-peak nuclei at densities < 10^^ 
g cm~^ [37|. Subsequently, when compressed to higher density due to accretion, these nu- 
clei undergo double electron captures ([A,Z] + 2^ ^ [A, Z — 2]) as /i^ increases, dripping 
of neutrons at p > pdrip ^ 6 — 7 x 10^^ g cm~^, and, when Z has been sufficiently reduced 
by the electron captures, pycno-nuclear fusions ([A,Z] + [A,Z] ^ [2A, 2Z], immediately fol- 
lowed by [2A, 2Z] + 2^ ^ [2A, 2Z — 2]). This sequence of nuclear reactions injects into the 
crust, depending on the details of the model, between 1.5-2 MeV for each accreted nucleon 
during its journey from the outer crust to the core. Known as "deep crustal heating" f38l, it 



constitutes the Qh heating term in Eq. 37 



In the Quasi-Persistent Sources the accretion outburst last long enough to heat the crust 
out of its initial thermal equilibrium with the core. However, since the core conductivity 
and heat capacity are very large its temperature does not vary appreciably during the 
outburst. Once accretion stops, the thermal relaxation of the crust can be observed, and is 
being continuously observed in the four known cases: 

- KS 1731-260: outburst of 12.5 yrs, with an estimated average accretion rate 
(M) - O.lMEdd, which ended in January 2001 |39|; 

- MXB 1659-29: outburst of 2.5 yrs, (M) - O.lMEdd, which ended in September 2001 EQI; 

- XTE J1701-462: outburst of 1.6 yrs, (M) - Mndd, which ended in August 2007 EB; 

- EXO 0748-676: outburst of 24 yrs, (M) - O.OlMgdd, which ended in August 2008 ll42ll : 

where M^dd ^ 10~^ Mq yr~^ ^ 10^^ g s~^ is the Eddington rate. We briefly present below 
two case studies. 
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Figure 11. Models for the crust relaxation of MXB 1659-29. See text for description. The 
six data points in the left panel are from HOl . assuming a source distance of 8.5 kpc. 



5.2.1. Mapping the thermal conductivity: MXB1659-29. 

The crust relaxation of MXB 1659-29 has been studied in detail by Brown & Gumming in 
||43]| . and our results amply confirm their analysis. The accretion outburst was long enough 
that the crust could reach a steady state: this is very important since it implies that the 
initial T profile for the crust relaxation was independent ofCy, providing some relief from 



the Cy/K degeneracy in Xth, Eq. 39 



As was shown in pSl there is a one-to-one mapping between the cooling curve, Te{t), 
and the temperature profile of the crust, T{z) at the end of the outburst at time ^o- At time 
t — t^ after relaxation commences, the observed surface temperature Te is determined by the 
temperature T{zt-to) ^ depth Zt-to such that the thermal relaxation time from the surface 



to this depth is Tth ^ t- (This is the "/^-effect" in Eqj39j) The schematic in the grey shaded 



inset in the left panel of Fig. [TT] shows: phase "1" when Te is determined by the outer crust 
evolution; in "2" it is controlled by the evolution of matter at densities p ^ 10^^ — 10^^ g 
cm~^; in phase "3" the evolution is sensitive to the deep inner crust; and, finally, in phase 
"4", the crust has relaxed with the core and Te reflects the core temperature. Approximating 
Cy and K by power laws in T, the evolution is described by power laws, i.e., straight lines 
in a Te-hog{t — to) plot. 

The thermal conductivity of a pure crystalline crust turns out to be much too high to 
reproduce observed cooling, but good fits are obtained when K is reduced due to additional 
scattering by impurities. The cooling curves in the left panel of Fig. [TT] illustrate three 
cases with impurity parameters gimp = 2.5, 5, and 7.5, as well as a pure crystalline crust, 
Qimp = 0. A value of Qimp = 5 is favored in this set of results, but is dependent on the 
assumed crust thickness and accretion rate (see f43l for a complete study). A finer study 
||45]| with a density dependent Qimp reveals that the cooling curves are mostly sensitive to 
the value of Qi^p at p > 10^^ g cm-^ so that MXB 1659-29, and also KS 1731-260, are 
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likely constraining the transport properties of the pasta phase. 

In the right panels of Fig. [TT] we show samples of T profile in the crust. In the upper 
panels are displayed the profiles at the end of the 2.5 yr accretion phase for the four models 
with various values of Qimp as labeled. In the lower panel we show the time evolution of 
the T profile for the model with Qimp = 5 and the gap model Bl after the end of accretion: 
times / — ^0 in days are indicated on the curves. 

An issue that was only briefly discussed in Ref. [43J is the range of uncertainty in 
the specific heat of the inner crust. Either the neutron Cy due to uncertainty the density 
dependence of the ^ So superfluid gap or the ion C\ from uncertainty in the effective mass A* 
and shear modulus of the ions can alter the evolution timescale. We illustrate the uncertainty 



on Cy in the left panel of Fig. 1 1 in the case Qimp = 5 by displaying four cooling curve for 
the four neutron ^Sq Tc curves of Fig|4j One sees a divergence of the curves with gaps 
A1/A2 from the ones with B1/B2 at the beginning of phase 2: this phase is sensitive to the 
physics just above the neutron drip point and gaps As result in a larger Cj) in this region, 
resulting in a slower cooling. During the later phase 3 the cooling curves with the gaps I's 
diverge from the ones with gaps 2's: this phase is sensitive to the physics in the deepest 
inner crust where the gaps Al and Bl vanish and imply a much larger CJ) than the gaps A2 
and B2. (The models with Qimp 7^ 5 where all calculated with the gap Bl.) The uncertainty 
on C\ is illustrated in the models with Qimp = 7.5: increasing A*, and thus Cy, resulting in 
only moderately slower cooling rates. (The models with Qimp 7^ 7.5 where all calculated 
with A* =A.) 

These results show that the dominant piece of physics controlling the evolution of 
MXB 1659-29 is the thermal conductivity (and, of course, the crust thickness), and sim- 
ilar conclusions hold for KS 1731-260 [44]: crust relaxation of cold neutron stars after 
a long accretion outburst are perfect laboratories to study the transport properties of the 
crust. 

5.2.2. A peculiar case: XTE 1701-462. 

Crust relaxation after the short, but strong, outburst of XTE 1701-462 presents distinctive 
behavior not seen in source with long and less powerful outbursts. XTE 1701-462 is quite 
orthogonal to MXB 1659-29 and KS 1731-260 in that its outburst was shorter, M about ten 
times larger, and its crust about 3 times warmer. It thermal relaxation is fit by a simple 
exponential with a decay time of about 100 days, in sharp contradistinction to the 500 days 
in MXB1659-29 and KS 1731-260 

The system is still cooling and, hence, not enough information is available to draw 
definitive conclusions, but some preliminary results are encouraging [461 . Two scenarios 



that fit well the observations to date are shown in the right panel of Fig. [12] Model "A" 
having a cold (^ 10^ K) core and "B" with a hot (^ 10^ K) core are both able to describe 
the data reasonably well. We note that the three hot points at days 225, 298, and 593 and 
marked by a "?" are most likely due to some short phase of residual accretion unrelated to 
thermal relaxation, (see BTIl for additional details). 



The two left panels of Fig. 12 show the evolution of the T profiles of these two models 
during the 1.6 yr long accretion outburst. These profiles should be contrasted with the T 
profile of MXB 1659-29 in the upper right panel of Fig. [TTJ it would have taken much more 
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Figure 12. Models for the crust heating and relaxation of XTE 1701-462 See text for 
description. The data are from BTIl . 

time for XTE 1701-462 to reach a steady T profile. For comparison, the "lOy" profiles in 
Fig. [12] correspond to a 10 yrs outburst, instead of 1.6 yr, and resembles more closely the T 
profile of MXB 1659-29 at the end of its outburst. The cusps in the early, 1 day and 1 week, 
profiles mark the locations of the energy sources: the major releases occur in the inner crust, 
from the pycno-nuclear reactions, while below pdrip only electron captures are present. The 
low density region, below 10^^ g cm~^, heats up rapidly due to its small specific heat while 
the region above it, but below Pdrip evolves more slowly due to its larger Cy and the absence 
of strong energy sources. The resulting dip in the T profile below Pdrip is precisely what 
leads to a much more rapid cooling in the observed light curve, compared to MXB 1659-29 
and KS 1731-260 where the heating was slower and the outburst long enough for heat to 
flow into this dip and produce a much smoother T profile. 

Keeping this region below Pdrip cold does require a low thermal conductivity, even lower 
than what was found in MXB 1659-29 (but as we noted above, this last system constrains 
K in the high density pasta phase, but not much in the Pdrip region). From the observation at 
day 175 till the last observation at day 1158 we see a plateau with a slow cooling: keeping 
the crust hot for such a long time also requires a significantly lower conductivity around 
and below Pdrip. 

While the inner crust of MXB 1659-29 at the end of the accretion outburst was at T ^ 
0.7 X 10^ K, in the case of XTE 1701-462 T reaches -2x10^ K: in the former case we 
are in a situation where impurity scattering strongly affected K and the ions to Cy was small 
(see Fig.[9]and|5]). Thus, not only does the astrophysical setting, high M and short outburst, 
but also the higher T encountered in XTE 1701-462 make it a peculiar system where the 
evolution is nearly independent of the electron capacity and electron-impurity scattering in 
the inner crust. The fourth system, EXO 0748-676, will present still another variant: very 
long outburst, low accretion rate, and high T ll46l . 
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6. Unresolved issues and future directions 

Renewed interest in the neutron star crust has already led to many recent developments 
and several of these are discussed in this book. In what follows we highlight a few areas 
where we can anticipate further progress and indicate why these improvements are needed 
to interpret observations. 

1. Existence and extent of the pasta phase: Since a large fraction of the crust (by mass) 
can be in the pasta phase it is important to determine the surface energy in asymmet- 
ric systems, and the density dependence of the nuclear symmetry energy. These are 
the key inputs needed to calculate the energy of the large lower dimensional struc- 
tures encountered in the pasta phase. These structures are favored over spherical 
nuclei only if the surface energy costs are modest. Calculations of neutron drops and 
asymmetric matter using realistic nucleon-nucleon interactions and ab-initio meth- 
ods such as Quantum Monte Carlo have already provided new insights. For example, 
calculations of the density profiles of neutron drops indicate that Skyrme functional 
used in the neutron star context typically underestimate the surface energy ll47]| . An 
exploration of different density dependencies of the symmetry energy, and the larger 
isospin gradient energies suggests a much reduced volume of pasta in the crust HSil . 

2. Velocity of shear modes in the inner crust: As we emphasize in this article at low 
temperature the velocity of the shear modes of the inner crust is a key input for 
both thermal and transport properties. The evolution of the shear modulus in the 
inner crust, especially in the pasta phases can differ greatly from those employed cur- 
rently. In addition, effects due to entrainment, finite nuclear size and polarizability 
can all tend to lower the shear speed in the inner crust. A self-consistent study of 
these effects will be useful to exploring connections between observations of thermal 
relaxation in a accreting systems and quasi-periodic oscillations in magnetars, since 
the shear mode plays a crucial role in both phenomena. 

3. Transport properties of the pasta phase: At high temperature, T MeV, the static 
ll49l and dynamic ll50l structure factors of the pasta phase have been calculated using 
molecular dynamics. However, since the regime of interest in accreting neutron stars 
is r <C < 1 MeV, transport properties will depend on the calculation of 5(0),^) at 
low temperatures where quantum effects and superfluidity are important. The shear 
modulus in the pasta region is still poorly known and will affect the velocity of trans- 
verse modes which was shown to play a crucial role in ^ and ^4] In addition, other 
collective excitations unique to lower dimensional pasta structures can be important 
and have been studied recently in the hydrodynamic limit lISTIl . More work is needed 
to understand how (al dente) pasta jiggles in order to calculate the heat capacity and 
dynamic structure factor at low temperature. 

4. How and when are Umklapp processes suppressed ?: The electron band gap is sen- 
sitive to short-distance details of the background lattice potential because ^pe ^ 
(5 — 6) / a. The effect of the nuclear form factor, impurities, and crystal imperfec- 
tions need to be studied to predict the transition temperature between normal and 
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Umklapp processes in the inner crust. This is especially important when interpret- 
ing the differences observed in thermal relaxation of systems with different accretion 
rates as they sample significantly different temperatures in the inner crust. Since im- 
purity scattering becomes important when T < Tum, observationally inferred values 
of 2imp may be erroneous if Umklapp suppression is incorrectly treated. 

5. Evolution of gimp in the inner crust: In accreting neutron stars the reaction paths 
ways are still not full explored in the inner crust. Here electron capture, pycno- 
nuclear reactions and nucleon emission, transfer and capture reaction that determine 
the abundances are very sensitive to poorly known nuclear structure effects such as 
pairing and shell gaps in extreme neutron-rich nuclei. It remains to be seen if generic 
expectations for gimp derived from nuclear structure studies are compatible with the 
inference that 2imp ^ 5 from the analysis of thermal relaxation observed in MXB 
1659-29. 

6. Are nuclear excitations relevant for T 10^ — 10^ K ? : Nuclear excitations are typi- 
cally ignored in the crust for T ^ 10^^ K. However, low-lying pygmy resonances, 
related collective surface modes, and perhaps even single particle excitations in the 
extreme nuclei encountered may affect both the thermodynamics and transport prop- 
erties if they occur at low energy <C COp . 
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